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Measurement of photons via conversion pairs in ^/snn = 200 GeV 
Au+Au collisions with the PHENIX experiment at RHIC 
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Abstract. Thermal photons can provide information on the temperature of the new state of matter created 
at RHIC. In the pr region of 1-3 GeV/c thermal photons are expected to be the dominant direct photon 
source. Therefore, a possible excess compared to a pure decay photon signal due to a thermal photon 
contribution should be seen in the double ratio (7/7(7r ))Measurcd/(7/7(7r°))simuiatcd, if sufficient accuracy 
can be reached. We present a method to reconstruct direct photons by measuring e + e~-pairs from external 
photon conversions. 



1 Introduction 



Direct photons are produced during all stages of heavy ion 
collisions at the Relativistic Heavy Ion Collider (RHIC). 
Because they do not interact strongly, they escape the 
medium unaffected by final state interactions and provide 
a promising signature of the earliest and hottest stage of 
the quark-gluon plasma (QGP) Q]. 

On a microscopic level, the main sources of direct pho- 
tons from a QGP are quark-gluon Compton scattering 
(<19 ~~ > 7?); quark-antiquark annihilation (qq — > 75) and 
bremsstrahlung involving thermalized partons |2j. Direct 
photons are also produced in initial hard scattering pro- 
cesses which involve the same reactions but among the 
incoming particles. 

At RHIC energies thermal photons are predicted to 
be the dominant source of direct photons in a px window 
between 1-3 GeV/c Q]. 

Direct photons have been measured with PHENIX in 
Au + Au collisions at ^/snn = 200 GeV |3] . The inclu- 
sive photon spectra measured with the Electromagnetic 
Calorimeter (EMC) have been compared to the expected 
background from hadronic sources, based on the measured 
7r and 77 spectra and a cocktail of other hadronic decays 
(7/, Kg, u>), assuming mj- scaling. 

Fig-fflshows the double ratio of the measured invariant 
yield ratio to the background decay ratio as a function 
of pt for minimum bias and for five centrality classes. 
The measurement of direct photons production at high px 
scales with the number of binary collisions in agreement 
with NLO perturbative QCD predictions and therefore 
confirms medium effects as the origin of jet suppression. In 
the low pr region, where a thermal signature is expected, 
a significant measurement remains limited by systematic 



uncertainties due to the energy resolution and the photon 
identification with the EMC. 

In order to overcome such limitations, dielectron pairs 
offer some advantages because of the superior resolution 
of charged particles at low momenta and excellent identi- 
fication of conversion photons; while other methods 0] try 
to use low mass dielectron pairs from internal conversions, 
the method presented here uses real photon conversions in 
the beam pipe. 
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Fig. 1. Double ratio of the measured invariant yield 
ratio, (7/ 7r °) M d > t° the background decay ratio, 
, as a function of pr for minimum bias and 

Background 

for five centralities of Au + Au collisions at ysjvjv = 200 GeV. 
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2 Thermal photon analysis 

The excellent capabilities of the PHENIX detector to mea- 
sure electrons suggest to circumvent the limitations of the 
conventional direct photon measurement [3j at low photon 
energies by measuring photons via their conversion pairs. 
The momentum resolution (~ 1 %) of the charged track- 
ing devices proves superior to the energy resolution of the 
EMC (~ 10 %) in the p T region of interest (1-3 GeV/c). 

Two steps are used to identify e + e~-pafrs from photon 
conversions. First a single electron identification cut is ap- 
plied, which require signals from at least two phototubes 
in the Ring Imaging Cherenkov Detector (RICH) match- 
ing to a reconstructed charged track in the Drift Chamber 
(DC). No further electron identification cuts were applied 
since the pair cuts (see Sect. 12. If) to separate conversion 
photons from other e + e _ -pairs are more efficient and pow- 
erful enough to provide a very clean photon conversion 
sample. 

The extracted photon conversions are tagged with pho- 
tons reconstructed in the EMC to determine the contri- 
bution from 7T° — > 77 decays (see Sect. 12.2)) . 

All yields are measured as a function of px of the e + e~- 
pair, which makes a direct comparison of the inclusive 
photon y: 
possible: 



photon yield, iVi, ncl , and the tagged photon yield, iY3" tag 
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-pair is already reconstructed, cancels. 



7 inc ' (Pt) 

7 hadr (p T ) 



(Pt) 



N 



! (Pr) 



Data 



jyhad,( pT) 

Hpt) 



(5) 



Sim 



The only remaining factor is the reconstruction efficiency 
of the photon in the EMC, e 7 (px), which has been deter- 
mined with Monte Carlo simulations to be, independent 
of the pair-p T , 95.0 ± 1.0 (syst) %. 



2.1 Photon Conversions 

Since the PHENIX tracking algorithm assumes the track 
to originate from the collision vertex, off- vertex conversion 
pairs are reconstructed with an artificial opening angle 
which leads to an invariant mass that is proportional to 
the radius at which the conversion occurs. 

Therefore, photon conversions that occur in the beam 
pipe material (Be, 0.3 % radiation length) at a radius 
of 4 cm are reconstructed with an invariant mass of ~ 
20 MeV/c 2 . Fig. [5] shows an invariant mass spectrum of 
e + e~-pairs in the range 0-0.1 GeV/c 2 . The peak from 
photon conversions in the beam pipe at 20 MeV/c 2 can be 
clearly separated from Dalitz decays ir° — -> 7e + e~, which 
dominate the spectrum below 10 MeV/c 2 , and combina- 
torial background pairs, whose contribution increases to- 
ward higher invariant masses. 



The measured yield of inclusive photons depends on the 
reconstruction efficiency e e + e - and the PHENIX accep- 
tance a e + e - of the conversion e + e - -pair . The tagged 
photon yield depends in addition on the efficiency to re- 
construct the second photon in the EMC £j(pt) and on 
the conditional probability / to find it in the EMC accep- 
tance, given that the e + e~-pair has been reconstructed 
already. Here £ 7 (pt), as well as all other yields and cor- 
rection factors, are determined as a function of the pt of 
the e + e~-pafr. Therefore, in a ratio of inclusive conversion 
photons to conversion photons which have been tagged as 
7T° decay products, the e + e~-pair reconstruction efficiency 
and acceptance correction factor cancel. 

A ratio of the hadronic decay photon yield, N^ adr , and 

the tagged photon yield from it decays, N* tag , is cal- 
culated with simulations, for which again the acceptance 
correction cancels. 
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Fig. 2. Invariant mass of e + e~-pairs before (black) and after 
(red) applying cuts on the orientation of the e + e~-pair in the 
magnetic field. 



iV 7 hadr ( PT ) = a e+e - 7 hadr (p T ) (3) 
Nf«* ( PT ) = a e+e - f 7 -° ( PT ) (4) 

The comparison of the ratio in data and in simulations 
in a double ratio leads to an expression that is equivalent 
to the ratio of inclusive and decay photons as shown in 
Eq. JSJ. In here also the conditional probability / of find- 
ing the second photon in the PHENIX acceptance, once 



The photon conversion pairs can be distinguished from 
Dalitz decays and purely combinatorial pairs by cutting 
on the orientation of the e + e~-pair in the magnetic field. 
The magnetic field inside PHENIX is parallel to the beam 
axis. Therefore conversion pairs, which have no intrinsic 
opening angle, are bent only in the azimuthal direction 
(i.e. perpendicular to the direction of the magnetic field). 
In contrast, the opening angle of Dalitz pairs and com- 
binatorial pairs can have any opening angle with respect 
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Fig. 3. Invariant mass of e + e~-pairs after applying pair cuts in 
the pair-pT region 0.8-1.2 GeV/c. The black line is indicating 
the remaining contamination by combinatorial background. 



to the magnetic field. Furthermore, the azimuthal opening 
angle of conversion pairs will always have the same sign, 
if one considers the ordered difference: 



Atpo = tpo (e~ 



<Po (e + ) 



(6) 



Fig. |21 shows the invariant mass spectra of e + e -pairs 
before (black) and after (red) applying these pair cuts. The 
yield from integrating the mass region < 35 MeV/c 2 of the 
conversion peak is corrected for the remaining p? indepen- 
dent contamination of 15.0 ± 2.0 (syst) % from combina- 
torial e + e~ -pairs which has been determined with mixed 
events (see black line in Fig. EJ - This leads to iVi, ncl (px), 
the total yield of conversion photons as a function of px- 



2.2 Tagging of Decay Photons 

To reveal which of these conversion photons come from 
n ° ~ * 77 decays, the e + e~-pairs in the conversion peak 
are combined with photons which have been measured 
in the EMC and their invariant mass is calculated (see 
Fig. UJ. Conversion photons that are identified as decay 
products of 7T° can be tagged as tag . This signal has 
a large combinatorial background due to the high photon 
multiplicity in Au + Au collisions. 

The combinatorial background is reproduced with an 
event mixing method, which creates uncorrelated pairs of 
photons and e + e~-pairs from different events. The mixed 
event spectra is normalized to the same event invariant 
mass spectra well outside the ir° mass region and sub- 
tracted. The normalization factor a is calculated as: 



1 ( FGl 

2 \BG1 



FG2 
BG2 



(7) 



with FGl (BG1) indicating the yield in the foreground 
(background) in the mass region 0-100 MeV/c 2 and FG2 
(BG2) the region 170-250 MeV/c 2 . 

The normalization factor has a relative statistical error 
on the order of 0.2 % which depends only on the size of 
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Fig. 4. Invariant mass of 7e + e -triplets in same events 
(black) and normalized mixed events (red) for e + e~-pairs with 
0.8 < p T < 5.0 GeV/c. 
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Fig. 5. Invariant mass of 7e + e~-triplets after background sub- 
traction for e + e~-pairs with 0.8 < pr < 1-2 GeV/c. A fit with 
a Gaussian is drawn and the resulting parameters shown in the 
box in the upper right of the graph. 



the normalization region in the foreground: 



1 



FGl + FG2 



(8) 



As an example, the resulting tt signal for e + e -pairs with 
0.8 < p T < 1.2 GeV/c is shown in Fig.0 

Mean and a are determined by a fit of the background 
subtracted data with a Gaussian. For comparison the data 
are also fitted to the sum of a second order polynomial 
and a Gaussian, to take into account the possibility that 
the shape is not completely described by the mixed event 
spectrum. The difference in the resulting mean and a is 
negligible. The mean and a obtained by the fit are then 
used to integrate the data in a region ±1.5 a around the 
mean. The integration region has been chosen to optimize 
the signal to background ratio. 

The statistical error on the extracted ir° signal is given 

by: 



4 



Y,FG(i) + aj:BG>(i) + (^Y, 
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Fig. 6. Number of conversion photons, A^ ncl (black), and num- 
ber of tagged 7T° , A^ tag (red), as a function of pr ■ (The yields 
are not corrected for the bin width in pr since this cancels in 
the ratio.) 
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Fig. 7. Shown in black is N^ adr , the number of e + e~-pairs 
from hadronic sources (ir° and if) as a function of pr from a 
fast Monte Carlo Simulation. The number of e + e~-pairs from 



7r" with the requirement that both, e 
were inside the PHENIX acceptance, A; 



e -pair and photon, 
tag , is shown in red. 



With FG(i) and BG'(i) being the yields in bin i of invari- 
ant mass spectrum in same events and normalized mixed 
events, respectively, the summations are performed over 
the integration region. It is important to note that the 
last term in Eq. © , is the square of the sum over the nor- 
malized background, and therefore, depends on the inte- 
gration region and is not bin independent. The systematic 
errors on the peak extraction have not been evaluated yet. 

In Fig. Elthe n° yield, A^° tag , as a function of pair- 
Pt is compared to the px distribution of all conversion 
photons A" lcl . Due to the limited PHENIX acceptance, 
7e + e~ triplets can not be reconstructed for e + e~-pairs 
with p T < 0.8 GeV/c. 



2.3 Simulations 

The contribution of hadronic decays has been determined 
with a fast Monte Carlo simulation of ir° and r\ Dalitz de- 
cays. A parameterization of the 7r° spectrum measured 
by PHENIX jS] has been used as input. The r] distri- 
bution has been generated assuming mx scaling {jpx ~ * 

Jp T + — irt^o) of the 7T° spectral shape and a nor- 
malization at high pr to -q/ir = 0.45 ± 0.10, which is in 
agreement with PHENIX data [5IT7] . 

In the pt region of interest, e + e~-pairs from Dalitz 
decays and from photon conversions have the same shape, 
which allows to circumvent a full Monte Carlo simulation 
of photon conversions in PHENIX, which would be a cum- 
bersome process due to the low conversion probability of 
0.2 % in the beam pipe. 

In order to produce inclusive and tagged photon yields 
comparable to the data, all particles (e + , e~ and 7) have 
been filtered into the PHENIX acceptance. Fig. F7| shows 
jyhadr^ p T distribution of all e + e~-pairs from hadronic 

sources, and tag , the e + e _ -pairs which originated from 
7T° decays. 



The ratio in Eq. 10, 7 incl (p T ) /j hadl (p r ), would result 
from dividing the ratio of spectra shown in Fig. |SJ 

(N™\ PT )/Nf^(p T )) 

V ' ' / Data 

corrected for the photon reconstruction efficiency, e 7 (px), 
by the one from the spectra in Fig. 

(N^( PT )/Nf^(p T )) . 

\ / Sim 

An excess above 1 would be interpreted as a signal from 
direct photons. 

3 Conclusions 

While a result can not be presented yet, the method in- 
troduced in here seems promising to find a signature of 
thermal photons in the low pt region as an excess in the 
double ratio as shown in Eq. ©. A detailed understand- 
ing of the combinatorial background in the invariant mass 
spectra of 7e + e _ -triplets which is currently the major 
source of systematic uncertainties will reduce the system- 
atic errors significantly. 
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